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Abstract. The goal of this paper is to derive the detailed description of 
the Enumeration Based Search Algorithm from the high level description 
provided in [16], analyze the experimental results from our implementa- 
tion of the Enumeration Based Search Algorithm for finding a regular 
bi-partite graph of degree 3, and compare it with known results from the 
available literature. We show that the values of m for a given girth g 
for (m, 3) BTUs are within the known mathematical bounds for regular 
bi-partitite graphs from the available literature. 



1 Introduction 

The goal of this paper is to develop the detailed description of the Enumeration 
Based Search Algorithm from the high level description provided in [16] and 
analyze the implementation results of the Enumeration Based Search Algorithm 
for finding a regular bi-partite graph of degree 3, and compare it with known 
results from the available literature, (m, r) BTU is our notation for a regular 
bi-partite graph that has been introduced in [1]. The high level description of 
the Enumeration Based Search Algorithm for searching a girth maximum (m, r) 
BTU has been described in [16]. The theoretical background behind BTUs has 
been introduced and explained in detail in [1] and [2]. 

2 Girth Maximization as a Extremal Graph Theory 
question 

We consider the problem of searching for a girth maximum (m, r) BTU as a 
question in Extremal Graph Theory by raising two related questions. 

1. Given girth g and r G N , what is the minimum value of m such that a (m, r) 
BTU has girth g . 

2. Given m, r G N; m ^ r, what is the maximum attainable girth for a (m, r) 
BTU ? 



2.1 Definitions 

We review definitions from [1] and [2] 
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Definition 1. {m,r) BTU 

A {m,r) Balanced Tanner Unit (BTU) is a regular bi-partite graph that can be 
represented by a m x m square matrix with r non-zero elements in each of its 
rows and columns. Every (m,r) BTU has a bipartite graph representation and 
an equivalent matrix representation. 

Definition 2. Girth maximum {m, r) BTU 

A labeled (m, r) BTU A is girth maximum if there does not exist another labeled 
(m, r) BTU B with girth greater than that of A . 

Definition 3. <?(/3i, /32, ■ • ■ , /3r-i) where Pi e P2{m) for 1 <i <r-l 

,p2,---, Pr-i) refers to the family of all labeled (m, r) BTUs with compatible 
permutations pi^P2, ■ ■ • ,Pr & Sm',Pi ^ C{Pi,P2, ■ ■ ■ ,Pi-i) for 1 < i <r that occur 
in the same order on a complete m symmetric permutation tree , xi^i < X2,i < 
. . . < Xr,i where pj = {xj^iXj^2 • • • xj^m)', ^ < j < r , such that Pi-i is the partition 
between permutations pi-i and pi for all integer values of i given by 1 < i <r . 

Definition 4. Optimal partition parameters for girth maximum (m, r) BTU. 

132, . . . , /3r-i S P2{b* k^^^) refer to optimal partitions derived in [2] such that 
there exists a girth maximum (to, r) BTU in <?(/3i, /32, ■ • ■ , Pr-i) ,where j3i refers 
to X)j=i~' b*k'^ = b*k'^~^ for 1 <i < r—1, with k gN obtained as a solution to 
b*k'^~^ = TO such that b €N is minimized. Thus, Pi, 02, ■ ■ ■ ,Pr-i G P2{b*k^^^) 
are Y!}^! b*k = b*k'-\ J2f^^' b*k^ ^ b*k'-^ , . . ., b * k"-^ = b*k'-^ 
and b*k'^~^ = b* k'^~^ respectively. 

2.2 Search for girth mELximum (m, r) BTU 

Search for a girth maximum (to, r) BTU refers to search for an optimal labelled 
(to, r) BTU in a family of labelled BTUs that we refer to as ^{Pi,p2, Pr-i) 
where Pi S P2{m) for 1 < i < r — 1. 



3 Girth Mciximization as a Extremal Graph Theory 
question 

We consider the problem of searching for a girth maximum (m, r) BTU as a 
question in Extremal Graph Theory by raising two related questions. 

1. Given girth g and r G N , what is the minimum value of m such that a (to, r) 
BTU has girth g . 

2. Given to, r e N; to r, what is the maximum attainable girth for a (to, r) 
BTU ? 
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4 Maximum Attainable Girth 

4.1 Maximum Attainable Girth for a (m, r) BTU 

We denote the maximum Attainable Girth for a (m, r) BTU as a function Qmax '■ 
{NU{0}}2^NU{0}. 

Theorem 1. The maximum attainable girth of a {m,r) BTU satisfies the in- 
equality gmax{rn,r) <2*k where k € N is obtained by minimizing b gN such 
that m = b* k'^~^ for r > 3 . 

Proof. Prom the optimal partition result from [2] for a (m, r) BTU for r > 3, 

we obtain that the maximum possible lenght of the maximum known cycle is 
2* k, where the optimal paritions are (3i refers to X^^Z^"* 6 * fc* = b* k^~^ for 
1 < i < r — 1 and k is obtained by minimizing 6 € N such that m = b* k^~^ for 
r > 3. We now need to show that gmax cannot equal 2 * fc for r > 3. This follows 
because of micro-partition cycles defined in [2] and their combinations which do 
not permit Qmax to equal 2 * for r > 3. Hence, the result follows. 

5 High Level Description Of Enumeration Beised Seeirch 
from [16] 

5.1 Enumeration Based Search algorithm for girth maximum (m, r) 
BTU for r > 3 

Wc find &, fc e N such that b is the smallest integer satisfying m = b* k'^~^; 
for( i = 2;i < r;i++) { 

Pi = Cj; min{b * — > b* such that (j, b * 6 * k^~^ — j) are 
relatively prime; 
if(i == 2 ) 

Pi-l = Ib*k'-'^ j 

else { 

Rearrange the {b * k'^~^,i) BTU such that pi-i = Ibi,ki-i ; 

Find qi-2 € Sf,*k'-^ such that it maximizes girth of (6 * BTU is formed 

by pi , . . . , pi-2 e Sfe^fei-i ■,Px = k * qa;-,! < X < i - 2; 

if(i != r - 1) 

Scale permutations Py = k * Qy-,! < y < i; 

} 

} 

5.2 Enumeration Based Search algorithm for a girth maximum 
(m, 3) BTU where m = b*k^ 

We find &, fc G N such that b is the smallest integer satisfying m = b* k"^; 
for( i = 2;i< 3;i++) { 

Pi = Cj; min{b * A;'~^ — > b* such that [j, b * k^~^,b * k^~^ — j) are 
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relatively prime; 
if(i == 2 ) 

Pi-l = Ib*k*-^', 

else { 

Rearrange the (6 * k'^~^,i) BTU such that pi^i = /^^fci-i; 

Find qi G Sb*k such that a girth maximum (b^k"^ , 3) BTU is formed hj p\,p2,Pz S 

Sb^k^-i;pi =k*qi; 

} 

5.3 Reorganizing the (6 * fe*~^, i) BTU such that Pi-i = Jj„fei-i 

Without loss of generality, we apply suitable permutations on depth and per- 
mutations labels on the {b * k'^~^,i) BTU in order to obtain pi-i = I^-tf^i-i. Per- 
mutations on depth and permutations labels have been explained and defined 
in [1] and preserve isomorphism since they correspond to row permutations and 
column permutations on the matrix representation of the (& * k^~^ , i) BTU. 

6 Detailed Description Of Enumeration Based Seairch for 
a girth mcLximum (fc^, 3) BTU 

To find permutation aqi € Sk such that a girth maximum (fc^, 3) BTU is formed 

by {pi,..,P3} { 

We enumerate all permutations qi with node at depth 1 fixed, such that partition 

between qi and 92 = Ik is {k) G P^ik) ; 

for(each enumerated permutation §1 ) { 

We scale up qi by k and P2 = Ik^ ; 

Ps = Cj where {j, fc^, fc^-j) are relatively prime; 

We compute the girth of this {P, 3) BTU; 

} 

We choose permutation qi that gives us the best girth; 

7 Detailed Description Of Enumeration Based Search for 
a girth maximum (m, r) BTU where 6, G N such that 
b is the smallest integer satisfying m — b * k^~^ 

To find permutations {qi, . . .,9^-2} G Si,*k'-^ such that a girth maximum (6 * 
k'-^,i) BTU is formed by {pi, ..,Pi} { 

We enumerate all permutations gi_2 with node at depth 1 fixed, such that par- 
tition between (7^-2 and q^^i = Ii^k^-'^ is (h * ; 
for (each enumerated permutation gj_2 ) { 

We permute {gi, . . ., qi-3} such that all partitions between any two permutations 

in the set {qi, . . ., qi-2, Qi-i} arc preserved; 

We scale up {qi, ■ . ., qi-2} by k and pi-i = Ib^.ki-i;Pi = Cj where (j, b * 6 * 
k^~^-j) are relatively prime; 
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We compute the girth of this (6 * A:* BTU; 
} 

We choose permutation qi-2 that gives us the best girth; 

8 Algorithm to Find Permutations Of {qi, . . ., qi-2} 

We permute {gi, . . ., qi-2} such that all partitions between any two permutations 

in the set {qi,. ■ ., qi-2, , Itrk^-^ } are preserved 

for( j = 2;i<&*fc*-2;j++) { 

d = Label at depth j of qi-2', 

Permutations On Depth {d,j); 

Permutations On Labels {d, j); 

We calculate the partition between permutations k * qi-2 and pi and girth; 

We accept the change to {gi, . . ., qi-2} if it improves the girth; 

} 

qi-i returns to Ib^:ki-2 after each run of the loop. 

9 Experimental Results for Implementation Of 
Enumeration Based Search 

Girth obtained for various values of m and for r = 3 has been shown Table 1 . 
We find that the values of m for a given value of girth g lie between the lower 
bound for m and improved lower bound for m from [13]. The execution time is 
too long for A; > 10 due to the algorithm being in EXPTIME. 

Table 1. Girth obtained for various of m and for r = 3 from Implementation 



k 


m 


r g 


5 


25 


3 8 


6 


36 


3 8 


7 


49 


3 10 


8 


64 


3 10 


9 


81 


3 10 


10 100 3 10 



10 Bound from [12] 

For q being a power of a prime fc > 3 , Lazebnik in [12] describes explicit 
construction of a g -regular bipartite graph on v = 2 * q'^ vertices with girth 
g>k + 5. 

If we consider this as a (m, r) BTU, we get r a power of a prime and m = 
r*^; fc > 3, girth g > log^(m) + 5. For g > 12, we obtain log^(m) > 7 which gives 
us m > r'' and we hence obtain m>3^ = 343 * 9 = 3087. 
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11 Lower bounds from [9] 

We quote the main theoom from [9], "Let G — (Vl, Vr, E) be a bi-partite graph 
of girth g — 2 * r, with til — {Vl)\ and nji — (Vr)|, the number of vertices 
on the left and right sides, and m = {E)\ the number of edges. Assume further 
that all vertex degrees in G are > 2 Then: hl > YA^oi^R^^'^'^'^'^K^L)^'""'^'^'^^ 
and UR > Elro'(^L)'="'^^/'H^K)^°°"^*/'^ where Ar = n„ey«(4-l)'^''/'", Al = 
TlveVL^'^^~^)'^''^"^ degree of vertex w." 

11.1 From another form of the bound in [9] 

From another form of the bound in [9], weobtainrii > '£lll{dR-lY'''^'^'/'^\dL-l)fl'''"'^'^'^^ 
and hr > Y^ZlidL l)='=^'(^/2)(rf^_i)i?oor(i/2)_por a (m,r) BTU with girth g, we 
obtain m > ^9/2-i(r-l)'='=»'(V2)(j._i)/ioor(»/2)_ 

Therefore, m > ^9/2-i(r-l)c«'KV2)+/ioor(z/2) 

For even integers i, ceil{i/2) + floor{i/2) = i 

For odd integers i , ceil{i/2) +floor{i/2) = (i + l)/2 + (« - l)/2 = i 
Therefore, m > Ef£o~'(^-l)' = - 1^'^-^+^ = ^"'(^^2)'^ 
Putting r = 3 and 5 = 12 we get m > ^ = 63. Putting r = 3 and g = 10 
we get m > ^ = 31. Putting r = 3 and g = 8 we get m > '"^^J^^ = 15 . 

11.2 Prom Main Theorem in [9] 

Derived from the main theorem. From [9], ul > Elro^(^fl)'"'''^^^(^i)-''°°''^'/^^ 
and riR > J21=o{AlT'''^^'^^H^r)^°'"'^'^^^ • For a (m,r) BTU with girth g , we 
obtain, Ar = {(r-l)'-/('"*'^)}™ = r-1 and Al = {(r-l)'-/^™*'-)}'" = r-1 . Thus, 
m > X;f£o"^(r-l)^'=*'(*/2)(r-l)^°°'-(^/2)_ Therefore, m > ^^'/2-l(^_l)ce»z(^/2)+/^oor(^/2) 

. For even integers i , ceil{i/2) + floor{i/2) = i . For odd integers i , ceil{i/2) + 
floor{i/2) = {i + l)/2 + {i - l)/2 = i . Therefore, m > Ef=o"^(^-l)' = 
(r - 1)9/2-1+1 ^ ^''"^^^^^"^ Putting r = 3 and 5 = 12 we get m > = 63 . 

Putting r = 3 and g = 10 we get m > ^ = 31 . Putting r = 3 and g( = 8 we 
get m > ^ 3^_~^ = 15 . 

12 Other Related Reseeirch 

Irregular LDPC codes with girth 20 in [11] and Regular LDPC codes of girth at 
least 10 from [10] . 
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13 Results from [15] 

Wc quote Theorem from [15] for even values of g since our current interest is 
only in bi-partite graphs. "For g >3 and 5 >3 put no{g, S) = '^*^^^ (^g__2) — 9 
is even. Then a graph G with minimal degree 6 and girth g has at least no{g, S) 
vertices." We use this result to compute no{g, S) for 5 = 3 and various values of 
g in Table 2 by simplifying the equation as no{g, 3) = 2 * {(2)^/-^-l} 

Table 2. Minimum value of no (5, 3) for different girths 5 for 5 = 3 from [15] 



g no{g,3) 

4 6 

6 14 

8 30 

10 62 

12 126 

14 254 



14 Results from [13] 

We quote theorems from [13] . 

1. "Given 6 > 3 and g > 3 , there exists a G", n < (2*(5)f with minimal degree 
of at least S and girth of at least g" . 

2. "Lower Bound n{g, 6) > if g is odd. n{g, S) > 

if g is even. Equality holds for 5 = 3 and g = {3,4,5,6,7,8} and g = 4:, 
(5 > 3" . 

3. "If g is odd, n{g+l,S) <2* n{g, S)" . 

4. "Upper Bound n{g,S) < ^^^^^^g^'"^^ ^ "(^'^) - ^*^^^~{5-2)'^~^^ 
if g is even" . 

5. "Let m > Y^fl^giS — ly = ^^~^l^_2) be an integer. Then there exists a S 
-regular graph of order 2 * m and girth of at least g" . 

6. "Most significant improvement of the bound for 5 = 3 , n{g,3) < 2^"^". 

15 Bound derived from [13] 

We derive the following bound from [13] , ^^^]j''_2 ^ ''^{Si ^) ^ '^*^''^s-2 — ~ 
for the minimum order n{g, 6) where g is its girth and 6 is its degree. By putting 
^ = 3, we obtain a simplified form of the above equation, (2)^/^ — 1 < n{g, 3) < 
4 * (2)^~^ — 1 which could be further simplified as 

29/2 _ 1 < ^^g^ 3) < 23 — 1. We calculate the bounds for (5 = 3 and the improved 
upper bound corresponds to n{g,3) < 2^~^ from [13] in Table 3. 
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Table 3. Lower Bound, Upper Bound and Improved Upper Bound for n{g,3) for 
different girths 5 for 5 = 3 from [13] 



9 


Lower Bound 


n{g, 3) Upper Bound 


n{g, 3) Improved upper bound n{g, 3) 


4 


3 


15 


8 


6 


7 


63 


32 


8 


15 


255 


128 


10 


31 


1023 


512 


12 


63 


4095 


2048 


14 


127 


16383 


8192 



16 Analysis for [17] and [18] 

From [17], we quote the following result, "If the degree is £> > 3 and girth 
g = 2 * r + l;r > 2, a simple lower bound for number of vertices of a regular 
graph is given by Uoig, D) = 1 + -^{{D-iy - I)." 

For D = 3 we simplify the equation as follows no(fif,3) = I + 3((2)(9-i)/2 _ 
While the exponent is similar to the lower bound in [13], we cannot apply the 
result as the girths take odd values and do not directly apply for bi-partite 
graphs. 

17 Analysis for [19] 

Wc analyze the girths obtained for various size of the matrices from [19] in Table 
4. However, these matrices have irregular degrees and hence a direct comparison 
with our obtained results might not be possible. 

Table 4. Girth obtained for various size of the matrices in [19] 



Girth Minimum N 

6 5 

8 9 

10 39 

12 97 



18 Analysis for [14] 

We quote from [14], "Ramanujan graphs X^''^ arc p+1 regular Cayley graphs of 
the group PSL{2,Z/qZ) if the Legendre symbol (2) = 1 and of PGL{2,Z/qZ) 
if the Legendre symbol (|) = —1 . X^''^ is bi-partite of order n = \{XP''^)\ = 
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q * {q^ — 1) and a bound on the girth is given by the equation, g{X'P''^) > 
41ogp(<z)-logp(4)". 

Putting p = 2 in order to get degree k — p -\- 1 — we obtain the inequahty 
g > 4:\og2{q) — log2(4) which can be simphfied as (g + 2)/4 > log2(g) in order 
to obtain 2'^3+2)/'i > ^ 

For each value of girth g, we calculate the minimum value of q such that q > 
2(s+2)/4 g^y^^^ ^j^g Legendre symbol (p/q) = —1 and then calculate n — q*{q^ — I) 
for p ~ 2 and degree fc = 3 in Table 5. 



Table 5. Analysis for [14] 



Girth 


min q,q> 2^3+2)/4^ f^^^^^ ^ 


— ln = q* (g^ 


— 1) Chosen p Degree k — p + 1 


6 


5 


120 


2 3 


8 


11 


1320 


2 3 


10 


11 


1320 


2 3 


12 


13 


2184 


2 3 



19 Conclusion 

Our implementation for the Enumeration based Search for a girth maximum 
(to, r) BTU finds the maximum attainable girth of a (m, r) BTU for r = 3 
and various values of to. The values of to for a given girth g are within the 
known mathematical bounds for regular bi-partitite graphs from the available 
literature. When we compare our results with bounds for more general graphs, 
or graphs with irregular graphs, a direct comparison may not possible since it 
is well known that for a given g and average degree, a lower number of vertices 
can be reached for irregular graphs. 
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